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Comparison of the ENEAR Peculiar Velocities with the 
PSCz Gravity Field 
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ABSTRACT 

We present a comparison between the peculiar velocity field measured from the EN- 
EAR all-sky D n — a catalog and that derived from the galaxy distribution of the 
IRAS PSCz redshift survey. The analysis is based on a modal expansion of these data 
in redshift space by means of spherical harmonics and Bessel functions. The effective 
smoothing scale of the expansion is almost linear with redshift reaching 1500km s _1 
at 3000 km s _1 . The general flow patterns in the filtered ENEAR and PSCz velocity 
fields agree well within 6000km s , assuming a linear biasing relation between the 
mass and the PSCz galaxies. The comparison allows us to determine the parameter 
j3 = f2 ' 6 /b, where f2 is the cosmological density parameter and b is the linear biasing 
factor. A likelihood analysis of the ENEAR and PSCz modes yields = 0.5 ± 0.1, in 
good agreement with values obtained from Tully-Fisher surveys. 
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1 INTRODUCTION 

In the standard picture for the formation of cosmic struc- 
tures via gravitational instability the peculiar velocity of a 
galaxy is generated by fluctuations in the mass distribution. 
For galaxies outside virialized systems, linear perturbation 
theory predicts 
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where Q is the mass density parameter, H a is the Hub- 
ble constant and 6 m is the mass density fluctuation field. 
If the relationship between the galaxy distribution, S g , and 
5m is approximately linear, S g = b5 m , then the parameter 
(3 — fi°' 6 /fo can be derived from the comparison between the 
observed peculiar velocity field and that predicted from the 
galaxy distribution. A particularly useful method for per- 
forming a velocity-velocity comparison is the modal expan- 
sion method developed by Nusser & Davis (1995, hereafter 
ND95). This method expands the velocity fields by means of 
smooth functions defined in redshift space, thus alleviating 
the Malmquist biases inherent in real space analysis. Fur- 



thermore, the modal expansion filters the observed and pre- 
dicted velocities in the same way, so that the smoothed fields 
can be compared directly. Because the number of modes is 
substantially smaller than the number of data points, the 
method also provides the means of estimating /3 from a like- 
lihood analysis carried out on a mode-by-mode basis, instead 
of galaxy-by-galaxy. The similar smoothing and the mode- 
by-mode comparison substantially simplify the error analy- 
sis. The modal expansion method has previously been used 
in comparisons between the 1.2 Jy IR A S predicted velocities 
and observed velocities inferred from Tully-Fisher (TF) mea- 
surements (Davis, Nusser & Willick 1996, hereafter DNW, 
da Costa et. al. 1998). In this paper, we perform a similar 
analysis using the recently completed redshift-distance sur- 
vey of early-type galaxies (hereafter ENEAR, da Costa et al. 
2000) and the IRAS PSCz redshift survey (Saunders et. al. 
2000). Because of differences in the nature of the data sets 
considered some slight changes in the method are required 
and are described below. Our goal is to investigate how well 
the velocity field mapped by early-type galaxies matches the 
velocity field inferred from the PSCz survey, and to obtain 
the parameter /3 yielding the best match. 
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In section 2, we briefly describe the ENEAR redshift- 
distance catalog. In section 3, we describe the modal expan- 
sion method as used here, present maps of the ENEAR and 
PSCz radial peculiar velocity field and perform a likelihood 
analysis to derive /3. A brief summary of our conclusions is 
presented in section 5. 



2 DATA 

We use a sub-sample extracted from the all-sky ENEAR 
redshift-distance survey (da Costa et al. 2000) comprising 
578 objects within cz < 6000 km s _1 , 355 field galaxies 
and 223 groups/clusters. Galaxies have been objectively as- 
signed to groups and clusters using redshifts taken from 
complete redshift surveys sampling the same volume. In- 
dividual galaxy distances were estimated from an inverse 
D n — a template relation derived by combining cluster data 
(e.g., Bernardi et al. 2000). The cluster sample consists of 
569 galaxies in 28 clusters. Over 80% of the galaxies in the 
magnitude-limited sample and roughly 60% of the cluster 
galaxies have new spectroscopic and photometric data ob- 
tained by the ENEAR survey. Multiple observations using 
different telescope/instrument configurations ensure the ho- 
mogeneity of the data. Furthermore, the sample complete- 
ness is uniform across the sky. 



3 THE MODAL EXPANSION 

An unbiased estimate of j3 = O 0,6 /b can be obtained from 
the comparison between smooth velocity fields with similar 
spatial resolution, derived from the ENEAR and PSCz data. 
To generate smooth fields we expand the peculiar velocities 
of both data in terms of smooth base functions. The expan- 
sion carried out here shares the general properties of that 
used by ND95, but differs in details. In their application to 
TF catalogs, ND95 defined Pi = 5 — log(I — m/si), where 
Si = czi is the galaxy redshift in km s _1 and Ui its radial 
peculiar velocity. The function P was then expressed by an 
expansion involving smooth functions. The final estimate of 
the smoothed velocity field was that obtained by minimiz- 
ing the scatter of the rotational speeds given the magnitudes 
in the inverse TF relation. The scatter was also simultane- 
ously minimized with respect to the the parameter of the TF 
relation. This led to an unbiased calibration of the inverse 
TF relation because the sample was mainly magnitude se- 
lected. The galaxy angular size and velocity dispersion in the 
D n — a relation do not uniquely fix the magnitude accord- 
ing to which the ENEAR sample is selected. So simultaneous 
minimization might lead to a biased estimate of the parame- 
ters of the D n — a relation. Although the bias is mild we use 
the calibration of the inverse D n — a given by Bernardi et al. 
(2000) by a regression of a on D n in clusters. We also express 
the peculiar velocity, u, rather than the function P in terms 
of smooth functions. Another difference is that ND95 used 
TF catalogs with all galaxies having the same relative dis- 
tance error which allowed an additional simplification in the 
application of the modal expansion method, namely, the ex- 
pansion in terms of orthogonal smooth functions. This made 
the TF velocity error covariance matrices diagonal. In the 
ENEAR sample, the relative distance error is not the same 



for all objects (galaxies and groups/clusters), so using or- 
thogonal functions does not offer any further simplification 
since the ENEAR error matrix remains non-diagonal. The 
lack of orthogonality slightly complicates the error analysis 
but does not affect the efficiency of the expansion. Choosing 
the spherical harmonics and Bessel functions to be our base 
smooth functions we write the radial peculiar velocity model 
as 

U(s, 6, (j>) = Qnim (knV(s)) - Cg] Yj m (6, <j>) . (2) 
l,m,n 

where the sum is over m = — I to +1, I = to l m am> an d 
n — to n m ax. For the reasons given in DNW, we formu- 
late our model to describe the velocity field with respect to 
the motion of the Local Group. The constant cn is non-zero 
only for the dipole term ensuring that u = at the origin. 
The function y(s) in the argument of the Bessel functions 
makes their oscillations match the radial distribution of the 
ENEAR data. Here we take y 2 = ln[l + ( j^,) 2 ], but other 
similar forms can be used as well. Also the expansion does 
not include a Hubble-like flow (u cx s) so we assume that 
any such flow has been consistently removed from the EN- 
EAR and PSCz velocities. The coefficients a„i m are found 
by minimizing 

X = ^2^ 2 [ui ~ u°f (3) 

where u° are the raw observed velocities and cn is the error 
of the velocity estimate resulting from observational uncer- 
tainties and intrinsic scatter in the D„ ~ a. For field galax- 
ies Oi = 0.23si and for groups of galaxies it is reduced by 
1 / v/ N g , where N g is the number of galaxies in the group. 



4 SMOOTH VELOCITY MAPS AND THE 
DETERMINATION OF j3 

We apply the modal expansion method to smooth the raw 
measured velocities of the 578 ENEAR objects within a red- 
shift of 6000 km s _1 (Bernardi et. al. 2000V We use 51 modes 
corresponding to l m ax = 4, Umax — 3 in (□). The smoothing 
scale of these functions is linear with redshift and matches 
the low resolution filter used in da Costa et al. (1998) (see 
their Figure 1). The smoothed velocities were then derived 
by minimizing (3) with respect to «; m „ assuming an error 
of Oi — 0.23si/-y/Ag in the raw velocities of the ENEAR 
objects. The reduced \ 2 P er d.o.f of the fit was 1.017, a sat- 
isfactory value in this type of analysis (see DNW, da Costa 
et al. 1998). 

Given an assumed value for f3 we interpolate the PSCz 
predicted velocity field, computed by Branchini et. al. 
(2000), to the positions of the ENEAR galaxies. Branchini 
et. al. obtained the PSCz velocities from the PSCz galaxy 
distribution with a Top-Hat window of width equal to half 
the mean particle separation at a given redshift. The PSCz 
fields are then expanded in the same orthogonal set of basis 
functions as employed for the ENEAR velocities. The PSCz 
and ENEAR velocities are guaranteed to have the same res- 
olution because the original smoothing of the PSCz density 
field is small compared to the resolution of the modal ex- 
pansion. 

The smoothed ENEAR velocities are shown in Figure Fl 
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in redshift shells 2000km s -1 thick. Comparison of this fig- 
ure and Figure 3 of da Costa et al. (1998) shows that the 
general flow pattern is remarkably similar. In the case of 
ENEAR, in the innermost shell very few prominent struc- 
tures are probed by bright ellipticals. However, in the next 
two shells a strong dipole pattern can be easily recognized, 
being of comparable amplitude to that of observed with the 
SFI galaxies. This dipole corresponds to the reflex motion 
of the Local Group, with infalling galaxies in the Hydra- 
Centaurus direction and an outflow towards the Perseus- 
Pisces complex. The quality of the match can be evaluated 
from Figure |i] which shows the residual velocity field ob- 
tained subtracting the smoothed PSCz field from that of 
the ENEAR, assuming f3 — 0.5. As can be seen the overall 
agreement is good with only a few more distant galaxies giv- 
ing large residuals. Note, however, that even though with a 
larger amplitude, the mismatch seen in the outermost red- 
shift shell at / ~ 0°, -60° ;Cb <; - 15° between ENEAR and 
PSCz correspond to mismatches in the comparison between 
SFI and 1.2 Jy IRAS velocity fields. This may correspond to 
a real mismatch between measured and predicted velocities 
which deserves further investigation. 

The filtered ENEAR and PSCz velocity fields are fully 
described by the modal expansion coefficients, a en and a ps , 
of the ENEAR and PSCz fields, respectively. Since the num- 
ber of these coefficients is significantly smaller than the num- 
ber of galaxies, it is more efficient to estimate /3 by compar- 
ing the modes rather than the individual galaxy velocities. 
As in da Costa et al. (1998) we define our best estimate 
for P as the value that corresponds to the minimum of the 
pseudo-x 2 



x 2 (/3)=y> 



4M] [T + M(/3)]- 



■(4) 



< Sa ] en Sal n > and M =< Sa ps Sa J ps > are the 



where T : 

the error covariance matrices of the coefficients a 3 en and a pg , 
respectively. For brevity of notation we have replaced the 
triplet n,l,m with one index j. The PSCz covariance ma- 
trix M incorporates errors due to (i) the uncertainty in the 
LG motion, which creates a dipole discrepancy between the 
ENEAR and the PSCz velocities, (ii) the discreteness in 
distribution of galaxies which propagates into the velocity 
field, and (Hi) small scale coherent (as in triple valued zones) 
nonlinear velocities that are not included in the PSCz recov- 
ered velocities. Details of how these error contributions are 
computed are in da Costa et al. (1998). Since the expansion 
functions are not orthogonal the ENEAR covariance matrix 
T has nonzero off-diagonal elements. This matrix is simply 



the inverse of 



where the derivatives are computed 



at the minimum of x 2 given by 

Given the covariance matrices, we compute the curve 
of the reduced X 2 (/3) as a function of f3, which is shown in 
the top panel of Figure ^. The curve was computed with 
an error of 150 km s _1 in the estimation of the LG motion 
and 160 km s _1 for the amplitude of nonlinear error in the 
PSCz field (see da Costa et. al. 1998). This amplitude of the 
nonlinear error was chosen to make the x 2 per d.o.f equal to 
unity at the minimum. In their analysis of the SFI and 1.2 
Jy IRAS, da Costa et al. (1998) obtained a lower value of 
90 km s" 1 for the amplitude of this error. The difference can 
be attributed, as expected, to a better match between the 



SFI and IRAS velocities and the increased nonlinearities in 
the PSCz velocity at the positions of the ENEAR galaxies 
which preferentially reside in high density regions. 

The minimum value of the x 2 is attained at f3 — 0.5, 
with the 1-sigma error being less than 0.1. We note that 
this result is not sensitive to the exact values adopted for 
the error estimates Another statistic indicating the good- 
ness of the match between the fields for various /3 is the 
correlation function of the residual u en — u ps between the 
smoothed ENEAR and PSCz radial velocities. This is shown 
in the bottom panel of Figure ^| for j3 — 0.2, 0.5, and 0.9. 
The amplitude of the PSCz field is small for f3 = 0.2, so 
the correlation function for this /3 is close to the correlation 
function of u en alone, while the opposite is true for (3 — 0.9. 
On the other hand, for /3 = 0.5 the correlation of the resid- 
ual velocity field is significantly smaller, indicating a good 
match between the measured and predicted velocity fields. 



5 SUMMARY AND DISCUSSION 

Using the modal expansion method of ND95 and the re- 
cently completed ENEAR redshift-distance survey and the 
PSCz redshift survey we have carried out a comparison be- 
tween the observed peculiar velocity field and that predicted 
from the distribution of PSCz galaxies. We find that the cor- 
responding smoothed fields agree well and the best match 
is obtained with j3 = 0.5 ± 0.1. This value is intermediate 
to those derived using the Mark III and SFI catalogs both 
based primarily on spiral galaxies. It is also consistent with 
the results obtained by Borgani et al. (2000) using an in- 
dependent method based on modeling the velocity correla- 
tion function. Note, however, that the discrepancy between 
the values determined from these methods and those ob- 
tained from the power spectrum analysis (e.g., Zaroubi et 
al. 2000) and density-density comparisons (e.g., Sigad et al. 
1998) still persist. The good agreement between SFI and 1.2 
Jy IRAS and between ENEAR and PSCz implies that the 
SFI and ENEAR velocity fields are also in good agreement. 
This suggests that the velocity maps obtained from the new 
distance- redshift surveys are a fair representation of the un- 
derlying velocity field, as the general characteristics of the 
observed flow fields are independent of the type of galaxies 
and distance indicators used. The good agreement among 
the values of (3 obtained using Mark III, SFI, ENEAR, 1.2 
Jy and PSCz catalogs gives further support to low values of 
(3 and point toward low-density cosmologies. 
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Figure 1. The sky projection in galactic coordinates as seen in 
the LG frame of the filtered ENEAR velocity field. 
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Figure 2. The residual velocity field (ENEAR minus PSCz) for 
f3 = 0.5 
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Figure 3. Top panel: curve of reduced pseudo-x 2 versus /3 com- 
puted using equation (5). Bottom panel: the correlation function 
of the velocity residual field for (i = 0.2, 0.5 and 0.9. 
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